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Water, the ubiquitous solvent, is also prominent in forming liquid-solid interfaces with cat- 
alytically active surfaces, in particular with promoted oxides. We study the complex interface 
of a gold nanocatalyst, pinned by an F-center on titania support, and water. The ab initio 
simulations uncover the microscopic details of solvent-induced charge rearrangements at the metal 
particle. Water is found to stabilize charge states differently from the gas phase as a result of 
structure-specific charge transfer from/to the solvent, thus altering surface reactivity. The metal 
cluster is shown to feature both "cationic" and "anionic" solvation, depending on fluctuation and 
polarization effects in the liquid, which creates novel active sites. These observations open up an 
avenue toward "solvent engineering" in liquid-phase heterogeneous catalysis. 



Highly dispersed gold nanoparticles supported on ox- 
ides have been shown to catalyze a number of impor- 
tant reactions, including low-temperature CO oxidation 
and the water gas shift reaction [l|, Q. Reducible ox- 
ides, inparticular titania (Ti02), are ideal catalytic sup- 
ports [31. The size of the gold particles substantially af- 
fects the catalytic activity, suggesting the key importance 
of metal/support interactions on a nanometer scale 
Reactions, and in particular CO oxidation, are believed 
to occur at specific active Au sites at the Au/Ti02 inter- 
face Although much is known regarding Au/Ti02 
catalytic activity in the presence of a gas phase, the com- 
plexity increases steeply when solvent is included. 

Liquid-solid interfaces as such are relevant to many 
industrial applications of great significance, such as 
(photo-)catalysis, solar cells, gas sensors, or biocompati- 
ble devices. In heterogeneous catalysis, it has been shown 
that the presence of water increases the observed rate 
of CO oxidation [l^. The degree of rate enhance- 
ment depends on the type of support used. In particu- 
lar for the case of Au/Ti02 catalysts it has been shown 
that their activity at about 3000 ppm II2O is so high 
that full conversion of CO is reached [l^. Thus, the 
Au/Ti02 surface displays a pronounced catalytic activ- 
ity toward the water-gas-shift (WGS) reaction. This fun- 
damental reaction represents a key reaction to produce 
extra H2 fuel from steam reforming, which is reversible 
and exothermic, according to the following reaction: CO 
+ II2O ^ CO2 + H2. The so-called carboxyl and redox 
mechanisms have been proposed for the WGS reaction on 
metal/oxide surfaces; in the former mechanism the CO 
species reacts with terminal hydroxyl groups, whereas in 
the latter CO reacts with an O atom from OH disso- 
ciation or from the oxide support. In both mechanism 
proposed the starting point is a H2O molecule that is 
initially adsorbed on the metallic cluster with its oxygen 
atom being attached to a metal atom. 

Moreover, even large (as opposed to nanoscale) gold 
particles, which are usually catalytically inert, show con- 
siderable oxidation activity at aqueous conditions [HI, 



12|. Because of its high dielectric constant, water may 
actively participate in chemical reactions for instance 
by stabilizing ionic species, thus speeding up reactions 
at liquid-solid interfaces, which opens up novel avenues 
to improve the performance of traditional heterogeneous 
catalysis at the gas-solid interface. Recently, it has been 
shown that the selective oxidation of alcohols in aque- 
ous phase over (Ti02, C)-supported Au catalysts is fa- 
cilitated by high pH conditions (T^]. This water- based 
approach to heterogeneous catalysis offers a sustainable, 
environmentally benign, and cheap alternative to tradi- 
tional processes that rely on both toxic and expensive in- 
organic oxidants and organic solvents IJ] . Despite these 
promising experimental findings, much remains unknown 
at the molecular level about the impact of water as a sol- 
vent on the reactivity of these modern catalyst systems. 

Several ab initio molecular dynamics (AIMD) [l^l 
studies have been recently reported focussing on the fun- 
damentals of the water-titania interface. For a water film 
on TiO2(110) it has been shown that at least two distinct 
layers form: molecules in the first layer arc sluggish and 
bind strongly to fivefold coordinated Ti sites, whereas 
those in the second layer interact only weakly with the 
substrate and diffuse rapidly [l6l | , thus yielding a highly 
anisotropic interface |17| . Moreover, water was not seen 
to dissociate at the coverages examined. Upon study- 
ing the thermodynamics of de/protonation of the rutile- 
water interface using free energy perturbation methods, 
a value of 0.6 eV was found for the dissociation free en- 



ergy of bulk water on defect- free rutile [18| . These stud- 
ies thus show that the ideal Ti02 surface is rather inert 
with respect to solvation by water. However, little is 
known about the interaction between water and defec- 
tive or metal-supporting Ti02 interfaces, being relevant 
to catalysis and industrial applications. 

Here, we present large-scale AIMD simulations aimed 
to investigate solvent effects at a gold nanocluster pinned 
by an F~center on TiO2(110) being in contact with 
liquid water. The aqueous solution is found to in- 
duce pronounced charge transfer and localization at the 
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nanocatalyst-liquid interface and stabilizes structures 
different from the gas phase. The TiO2(110) surface has 
been modeled by four 0-Ti202-0 tri-layer (4x2) super- 
cell slabs separated by more than 15 A. The most com- 
mon point defects on the TiO2(110) rutile surface are 
oxygen vacancies in the twofold coordinated O rows Q . 
Therefore, a Aun nanocluster was grown on TiO2(110) 
where an F-center created by a surface O vacancy acted 
as anchoring site for initial Au nucleation; note that 
this blocks the O defect to interact with water and thus 
prevents splitting of water molecules at this F-centcr. 
The lowest energy structure of Aun adsorbed on a re- 
duced rutile surface and employed in the present study is 
shown in Fig. [TJc)). In order to create the Auii/Ti02- 
water interface, the space between the slabs has been 
fully filled with 53 II2O molecules (see Fig. [TJa)). It 
has been shown that Aur is the smallest Au cluster 
on rutile TiO2(110) with a measurable reaction rate for 
CO combustion [19(. The Aun model, although be- 
ing a simplification of the nanometersized clusters on 
TiO2(110), appropriately mimics the active sites located 
at the nanogold/oxide interface where the oxidation pro- 
cess takes place @, [2^ . 
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FIG. 1: Ball and stick model of the Aun nanocluster 
pinned by an F— center on the TiO2(110) oxide surface in 
contact with liquid water (a) and with a gas phase (c). 

Rod, big blue, small blue and brown spheres are O, Ti, H and 
Ti"^"*" sites, respectively, whereas the surface O vacancy is shown in 
green. Violet, black and yellow spheres correspond to top, middle 
and bottom Au sites of the Aun nanocluster, respectively. Time 
evolution of the Lowdin charges of the top (violet lines) and middle 
(black lines) Au atoms of the supported gold nanocluster in liquid 
water (b) and in the gas phase (d). Red and orange triangles corre- 
spond to reference Lowdin charges computed for a single Au+ (aq) 
cation (10.924 |e|) and Au°(aq) atom (11.210 |e|) in liquid water, 
respectively. 



The binding of the Aun cluster on the reduced 
TiO2(110) support entails a strong charge rearrangement 
at the metal/oxide contact (see Supporting Figure 1), 
the adsorption energy being —2.19 eV. In the case of an 
isolated O vacancy on the stoichiometric TiO2(110) sur- 
face, the charge neutrality is maintained by the presence 
of two Ti'^"'" ions thus creating an F~center. Spin density 
and bonding charge analyses reveal that the Aun clus- 
ter, once adsorbed, leaves a reduced substrate with three 
Ti^"*" ions, and ^ 0.4 |e| are transferred from the metal 
cluster to the oxide. This results into a slightly positively 
charged Au^f cluster supported by a reduced TiO2(110) 
oxide surface. 

In order to reveal the solvent-induced charge re- 
arrangement at the metal-liquid contact, we have 
performed PBE + U AIMD simulations of the 
Auii/TiO2(110) nanocatalyst in aqueous solution. 
Additionally, corresponding gas phase simulations have 
been carried out to provide the solvent-free reference 
situation. In both cases the finite temperature dynamics 
preserves three reduced Ti'^^ sites (see Supporting 
Figure 2). Snapshots of the AIMD simulations have 
been collected every ~ 0.2 ps thus generating a set of 
representative configurations for electronic structure 
analyses. In the following, we will refer to Au atoms of 
the supported Aun nanocluster as top, middle and bot- 
tom sites, which are depicted as violet, black and yellow 
spheres in Fig. [U respectively. A detailed investigation 
of the interaction between II2O molecules and the bare 
oxide surface far from the supported metal cluster is 
out of the focus of this study , but we note in passing 
that our AIMD simulations are in agreement with those 
of [iB-Q in the regions. In particular, water at the 
liquid-oxide contact is slow, strongly bonded to fivefold 
coordinated Ti sites, and is not seen to dissociate. 

The charge dynamics of the Au sites was extracted by 
computing, as a function of time, the Lowdin charges, 
see Fig.[Tl where violet and black lines correspond to the 
charge evolution of top and middle Au atoms, respec- 
tively. As demonstrated by Fig. [TJb) , two distinct charge 
patterns, corresponding to top and middle Au sites, can 
be clearly distinguished in aqueous solution; note that 
Au atoms in direct contact with Ti02 (yellow spheres 
in Fig.[TKa)) cannot be solvated by the liquid phase and 
thus do not exhibit a clear trend. The average charges 
on top and middle Au atoms are 11.028 ± 0.006 and 
11.220 ± 0.006 |e|, respectively (see Supporting Table 1 
for statistical analysis). It is noted in passing that a very 
similar charge separation has been found for a gold clus- 
ter of different size and shape using an Aui3/TiO2(110) 
nanocatalyst in contact with water. Interestingly, at 
around 6 ps a considerable amount of charge is trans- 
ferred to a specific top Au site due to solvent fluctua- 
tions, which reaches a value that is typical of middle Au 
atoms, see Fig. [IJb). In stark contrast, no such charge 
separation and fluctuation effects can be identified in the 
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FIG. 2: Dynamics of the water ^ nanocluster charge trans- 
fer (CT). Violet and black lines correspond to top and middle Au 
atoms, respectively 
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FIG. 3: Water-induced charge flow. Left and central 
panels: charge difference, Ap(r), for a representative snap- 
shot of the Auii/TiO2(110)-water interface. Electron de- 
pletion/accumulation is depicted by blue/red isosurfaces at 
±0.02 |e|/A^. Right panel: charge difference integrated in planes 
perpendicular to the surface, Ap{z), as a function of the height, z. 



gas phase reference dynamics in Fig. [TJd). We conclude 
that the hquid phase promotes and stabihzes a net charge 
separation in the supported Au cluster, where top and 
middle atoms can be clearly distinguished on the basis 
of their electronic structure. As a result, water steers 
the Au sites towards either one or the other preferred 
charge value, depending on their specific location within 
the nanocluster and going hand in hand with different 
solvation patterns, vide infra. 

To further characterize this puzzling process we have 
quantified the amount of charge transfer (CT) from 
Au atoms to the aqueous solution and back using Badcr 
analysis j " 
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For this purpose we have computed 
the difference between the Bader charges of the Au 
atoms in the fully solvated Auii/TiO2(110)-water sys- 
tem and those of the corresponding Auii/TiO2(110) gas 
phase system, where the latter is obtained upon remov- 
ing the solvent while keeping all atomic positions fixed. 
As depicted in Fig. [2] we find that top Au sites trans- 
fer charge to the solvent, with an average CT of about 
0.13 ± 0.05 |e|, while middle Au atoms attract about 
0.06 ± 0.03 |e| per Au atom from water. 

In view of these pronounced CT effects, we probe 
solvent-induced morphology changes of the metal clus- 
ter with respect to the gas phase. A set of trajectory 
configurations of the solvated nanocatalyst has been se- 
lected where water has been removed before quench- 
ing and optimizing the remaining Auii/TiO2(110) sys- 
tem. In order to generate the proper gas phase reference, 
Auii/TiO2(110) in the absence of water has been run 
at the same temperature of 450 K before applying the 
same quenching protocol; note that this is the relevant 
experimental temperature used for selective liquid-phase 
alcohol oxidation in aqueous solution using gold/titania 
catalysts. In case of the quenched solvated interface, 
the same local minima and no significant change in the 
Lowdin charges were observed as in the solvated state. 
However, the gas phase reference system yielded several 
local minima with significant changes of both the real- 



space structure and the Au charges (see Supporting Fig- 
ures 3 and 4) . These findings show that the liquid phase 
stabilizes structures and charge states within nanogold 
that arc distinctly different from those observed in the 
gas phase, which eventually induces novel morphologies 
and active sites, respectively, due to solvation. 

Let us now investigate the solvent-induced electronic 
charge redistribution at the nanocatalyst. In the fol- 
lowing, we refer the electronic charge density of the sol- 
vated Auii/TiO2(110)-water system to that of the iso- 
lated one, Auii/TiO2(110), and the separated solvation 
water. We thus compute the difference, Ap{r), between 
the charge densities of the combined system and the sep- 
arated components while freezing all atomic positions. 
A representative example of water-induced charge fiow 
is depicted in Fig [3l where the positive (negative) com- 
ponent of the charge difference Ap_|_(_)(r) is shown by 
red (blue) isosurfaces. This reveals that charge depletion 
(blue) occurs mainly on top Au atoms, whereas addi- 
tional electron charge (red) prefers to accumulate on mid- 
dle Au sites. Furthermore, charge depletion and accu- 
mulation visible at the oxide- water contact can be traced 
back to the interaction between water and the TiO2(110) 
oxide support. Here charge depletion occurs mainly on 
top of surface O atoms to which H2O molecules point 
with their H atoms, whereas charge accumulates close 
to five-fold coordinated Ti atoms, being in turn coor- 
dinated by water O atoms. This is quantified by con- 
sidering the charge redistribution perpendicular to the 
support, Ap{z) (see Fig. |3]). Here Ap{z) is the sum of 
the positive and negative charge density components and 
integrated in planes perpendicular to the surface. This 
quantitative analysis reveals that, on the one hand, there 
is a net charge transfer of ~ 0.22 |e| from top Au atoms 
to H2O molecules (obtained by integrating Ap{z) within 
the shaded orange area in Fig.|31), leading to slightly pos- 
itively charged Au''+ atoms at the corresponding sites. 
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FIG. 4: Snapshot of a typical hydration structure for 
the Auii/TiO2(110) nanocatalyst in aqueous solution. In 

the cationic/anionic solvation pattern H2O molecules (violet /black 
spheres) are pointing with their O/H atoms toward top/middle 
Au atoms (violet/black spheres) (a), (b) Time evolution of the 
Au-O distances of the gold atoms depicted in violet. 



On the other hand, there is a reverse charge flow of 
^ 0.37 |e| from water to middle Au atoms (see yellow 
area in Fig. [3]) , thus creating partially negatively charged 
Au''~ atoms upon solvation. Similar results were ob- 
tained for all the representative structures sampled along 
the trajectory. This can be compared to the average neg- 
ative and positive CT obtained from the Bader analysis. 
In total, this amounts to a structure-specific charge ex- 
change between metal and water of —0.34 and -fO.46 |e|, 
in qualitative agreement with the charge difference anal- 
ysis. As a result, top and middle Au sites experience 
opposite charge flow upon interaction with the solvent, 
as clearly shown in Fig. [3l which connects to a peculiar 
solvation pattern. 

We thus conclude by discussing the solvation pattern 
around the Auii/TiO2(110) nanocatalyst, see Fig. |4l[a), 
which is found to feature both so-called "cationic" and 
"anionic" solvation depending on fluctuation and polar- 
ization effects. This is similar to what we have observed 
previously for a single Au"'"(aq) cation and neutral atom, 
Au°(aq), in aqueous bulk solution [25| . In particular, the 
Au"*" aqua ion forms a quasi-linear molecular structure, 
II2O-AU+-OII2, in which two water molecules arc firmly 
bound through their oxygens to the Au^ site, as expected 
for a cation, yielding an Au+-0 distance of ~ 2.04 A. 
Upon adding an electron to Au+ to yield Au°, we found 
that its solvation pattern changes distinctly. The lat- 
ter features two distinct regions, supporting a Janus-type 
behavior. In one region a typical anionic solvation pat- 
tern prevails where two water molecules are pointing with 
their H atoms toward the metal. In other regions water 
binds via its O site to the Au°(aq), thus giving rise to a 
typical cationic solvation pattern. 

Around the nanocatalyst, two distinct solvation pat- 
terns corresponding to different Au charge states can be 
identified. The solvation around middle Au atoms con- 



sists of H2O molecules pointing with their H atoms to- 
ward the Au atoms like in an anionic solvation shell. 
In addition, several water exchange processes are ob- 
served as for Au°(aq) in the bulk [25[. At variance with 
those sites, the solvation around top Au atoms features 
a typical cationic solvation pattern, where a single H2O 
molecule point to a top Au atom via its O site. A pro- 
jected density of states analysis reveals that the interac- 
tion between gold atoms and solvating water molecules 
results from the overlap between the Au(d) and 0(p) 
states of the respective gold and oxygen sites. Interest- 
ingly, no H2O exchange processes are observed in this 
case again in line with the Au"*" (aq) cation case [l^l . As 
depicted in Fig. [U the average Lowdin charges of top 
and middle Au sites are in qualitative agreement to those 
computed for a single gold cation (10.924 |e|) and neutral 
gold atom (11.210 |e|) in aqueous bulk solution. Based 
on this comparison, we are in a position to explain the 
sudden switch of the Lowdin charge involving a specific 
top Au atom and taking place at ^ 6 ps in Figs. Wlh) 
and [2] In Fig. (H^b) we plot the distances between the 
top Au atoms with typical "cationic" solvation pattern 
and the O atom of its solvating II2O molecule. This 
demonstrates that two Au atoms are always "bonded" 
to the the same II2O molecule along the full trajectory, 
the time-averaged distances between Au and O sites be- 
ing 2.25 and 2.35 A. The situation is very different at 
the third top Au site. This Au atom is initially cation- 
ically solvated yielding a Au~0 distance around 2.3 A. 
At ^ 6 ps the H2O molecule rearranges and points to 
the same gold atom with one of its H atoms, a fiuctu- 
ation which induces the change of charge state of this 
Au site discussed earlier. This reveals the intimate con- 
nection between the creation of excess charge at specific 
nanoclustcr sites and the preferential orientation of the 
intcrfacial water molecules. As discussed in the introduc- 
tion for the carboxyl and redox mechanisms proposed for 
the WGS reaction on metal/oxide surfaces, the first step 
involves II2O molecules that are initially adsorbed on the 
metal cluster such that their O atoms are coordinated to 
an Au site. Therefore our simulations suggest that top 
site Au atoms, which feature a typical "cationic" solva- 
tion pattern at the water-gold interface, are the most 
promising candidates for being the active site for CO 
oxidation via the WGS reaction, whereas "anionically" 
coordinated Au sites arc expected to be chemically inert. 
The very same argument will hold for other liquid-phase 
catalytic reactions - thus providing a novel mechanism 
to create active sites at catalyst-water interfaces in much 
more general terms. 

In conclusion, finite temperature PBE + U simulations 
have been employed to investigate solvent-metal interac- 
tions on the gold/titania nanocatalyst in liquid water. 
Comparing our liquid phase to gas phase and bulk sol- 
vation reference data demonstrates that interfacial water 
alters both the structure and electronic properties of the 
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supported metal cluster in a significant way, including the 
creation of active sites. The fundamental phenomena ob- 
served will have a profound impact on understanding the 
role of solvent in heterogeneously catalyzed liquid-phase 
reactions, where de- and re-solvation processes play a key 
role. 

Methods 

All calculations have been performed using spin-polarized 
PBE + U [26ll , and ultrasoft pseudo pot entials [27l | as implemented 
in CPMD [28l and Quantum Espresso |29ll . In line with our previous 
work I21I. |22|. I30II . the value of U = 4.2 eV was adopted. We have 
carefully checked the charge distribution dependence on U by re- 
computing the Lowdin charges of the top and middle Au atoms of 
the supported gold nanocluster in liquid water for selected snap- 
shots along the trajectory shown in Fig. [TJb) using four different 
U values and found no significant change (see Supporting Material 
for data). 

The AIMD simulations [l5l | used Car-Parrinello propagation, 
the canonical ensemble was established at 450 K with a Nose- 
Hoover chain thermostat, and about ~ 10 ps trajectories have been 
generated for analysis (see Supporting Material for details). 
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